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ABSTRACT 

We present the formalism and numerical results for torsional oscillations of relativistic 
stars endowed with a strong dipole magnetic field, assumed to be confined to the crust. 
In our approach, we focus on axisymmetric modes and neglect higher-order couplings 
induced by the magnetic field. We do a systematic search of parameter space by 
\^ ' computing torsional mode frequencies for various values of the harmonic index t and 

^SJ , for various overtones, using an extended sample of models of compact stars, varying 

in mass, high-density equation of state and crust model. We show that torsional mode 
frequencies are sensitive to the crust model if the high-density equation of state is 
very stiff. In addition, torsional mode frequencies are drastically affected by a dipole 
, magnetic field, if the latter has a strength exceeding roughly 10^^ G and we find that 

the magnetic field effects are sensitive to the adopted crust model. Using our extended 
numerical results we derive empirical relations for the effect of the magnetic field on 
torsional modes as well as for the crust thickness. We compare our numerical results 
2 . to observed frequencies in SGRs and find that certain high-density EoS and mass 

4— » ' values are favored over others in the non-magnetized limit. On the other hand, if the 

, magnetic field is strong, then its effect has to be taken into account in attempts to 

formulate a theory of asteroseismology for magnetars. 
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1 INTRODUCTION 



There is growing evidence that Soft Gamma Repeaters (SGRs) could be magnetars experiencing starquakes that are connected 
(through the intense magnetic field) to gamma ra y flare activity. Magnetars ar e thought to be neutron stars with very strong 
magnetic fields, greatly exceeding ~ 5 x 10^^ G jDuncan fc Thompsonlll992h . The fiare activity consist of an initial short 
peak in the hard part of the spectrum, followed by a decaying softer part (tail) which last for hundreds of seconds. Recent 
analysis of detailed observations from several SGRs ( Israel et al. 20051 : Strohmaver fc Wattdboosl : Watts fc Strohmaver 20061 ) 
revealed that the decaying part of the spectrum exhibits a number of oscillations with frequencies in the range of a few 
tenths of Hz to a few hundred Hz. There are three events detected up to now which can be associated with crust oscillations 
of m agnetars. The first event was detected alr eady in 1979 from the source SGR 0526-66 jMazets et al.lll979l : iBarat et al 



1983ll , the second in 1998 from SGR 1900-1-14 (iHurlev et al.lll999l). while the third and most energetic one was observed in 
December 2004 from the source SGR 1802-20 (jTerasawa et al.1 20051 : IPalmer et al.ll2005h . Analysis of the tail os cillations of 



SGR 1806-20 revealed the presence of oscillations at appro ximately 18, 26, 29, 92.5, 150, 626.5, and 1837 Hz ijisrael et al 
20051 : IWatts fc Strohmaver ' 2006: Stroh maver fc Wattsl[20oi ). for SGR 1900-1-14 the detected frequencies were 28, 54, 84, and 
155 Hz~(lstrohmayer fc Watta.2005 '). while for the case of SGR 0526-66 there was only one frequency identified at 43.5 Hz 
I Barat et al.lll983l ). 
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Magnetar starquakes may be driven by the evolving intense magnetic field, which accumulates stress and eventually leads 
to crust fracturing. The excited seismic oscillations are known to be of two types: shear modes, which are polar- type oscillations 
and torsional modes, which are axial-type oscillations. The latter are thought to be more easily excited during a fracturing 
of the crust, since they only involve oscillations of the velocity, with the shear modulus of the crust acting as the necessary 
restoring force. Slowly rotating stars still maintain a nearly spherical shape and the velocity field of torsional oscillations is 
then divergence-free, with no radial component. Torsional modes (t-modes) are labeled as itn, where I is the angular index, 
while the index n corresponds to the number of radial nodes in the eigenf unctions of the overtone s for a specific I. Sh e ar an d 
torsional modes ha v e been studi e d ma i nly in New tonian theory, see e.g. [Hansen fc CiofB l|l980h : iMcDermott et al.1 dlQSSl ): 



Carroll et al 



( 198611 :IStrohmaverl (Il99lh 



i Schumaker fc Thorndll983l : 



Leini 



1994 



Le3 ll2006l). w hile there are only a few studies of tor sional modes in general relativity 
Messios et al. 2001 ; Samuelsson fc Andersson 20061 ) . 



In the Newtonian, non- magnetized limit Hansen fc CiofB 1 1980l ) found that the period ^Pq of the fundamental torsional 



modes depends mostly on the radius of the star, R, the speed of shear waves, Us, and the angular index i, via 

tPo^2nW+l)]~^'^R/us, (1) 

where Us = (/i/p)^''^, with /i and p being the shear modulus and the density, respectively. Furthermore, it was found that the 
period of overtones of index n is essentially independent of the angular index I (provided / does not greatly exceed n) and is 
basically determined by the crust thickness Ar: 



iPn « 2n ^Ar/us- 

From the above two relations it follows that the relative crust thickness, Ar/R, is given by 
Ar Tvn ifo 



y^£{£ + 1) if" 



(2) 



(3) 



which is independent of the details of the equation of stat e in the crust and basically only depends on the ratio of the frequency 
of the overtones to the fundamental frequency, see also ISamuelsson fc AnderssonI ((2006). It thus becomes obvious that the 
successful identification of both the fundamental torsional frequencies and their overtones could allow for the determination 
of the crust thickness, which in turn could lead to information on the high-density part of the equation of state. In such a 
determination, of course, general relativistic and magnetic field effects on the frequencies of torsional modes should properly 
be taken into account, which is one of the aims of the present paper. 

Relativistic effects have been found to significantly increase the fundamental £ — 2 torsional mode period by roughly 
30% for a typical M = IAMq, R = 10km model. rPuncan. (,1998i) derived the following empirical formula for the period of the 
fundamental £ = 2 torsional mode of a relativistic star 

-2 



fPo(ms) ^ 33.6P1 



0.87 + 0.13Afi.4Pr, 
(1.71 -0.71Mi.4Pro')'''^ 



(4) 



where -Rio = -R/(10km) and Mio = M/{1.4:Mq). In this relation, two effects have been taken into account, the (relativistic) 
gravitat i onal r edshift an d the (Newton i an) ac celeration of gravity g = GM/R? at the surface. More details can be found in 
Duncan ( 1998l l. s ee also Messios et al. ( 2001 '). In addition, the formula Q agrees well with the fully relativistic empirical 
relation found bv I Messios et al.l (|200ll ) 



<Po(ms) « 34Pi 



(5) 



The latter is an extension of the Newtonian Hansen fc Cioffi formula, based on a particular 1. 4M(7) model. On the othe r 
hand, the effect of a strong magnetic field, B, on torsional modes was considered by DuncanI ( 1998h and Messios et al. ( 2001 ). 



Assuming that the shear modulus /i is augmented by the magnetic field tension, _B^/47r, one can derive the following estimate 
for the period of a torsional mode 

P«p('"[l + (B/B^)]-i/\ (6) 

where P'*^' is the period in the limit of vanishing magnetic field and — ^f^Kji. However, Messios et al. ( 200 ll ) found that if 
the structure of the magnetic field is close to monopolar in the crust, the period of the fundamental £ = 2 torsional mode does 
not follow the above relation. This underlines the importance of properly taking into account the structure of the magnetic 
field, so that Equation (O must be generalized to particular magnetic field multipoles. 

An essential ingredient in computing torsional mode frequencies is the assumpt ion one makes about the shear modulus 
p. Here, we adopt the zero-temperature limit of the approximate formula derived bv lStrohmaver et al.l (|l991f ) 



0.1194- 



(7) 



where rii is the ion number density, c? — Z/A-KJii is the average ion spacing, and -fZe is the ion charge. Notice that this 
approximate formula for the shear modulus is derived on the assumption of a bcc crystalline state, which is averaged over 
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all directions. For the detailed composition o f the crust, we consider a m odern EOS bv lDouchin &: Haensell 1 200 ll ) (hereafter 
DH) and, for comparison, an older EOS by iNegele fc VautherinI (jlQTd 'l (hereafter NV). The two EOSs for the crust differ 
both in composition as well as in the value of the density at the base of the crust. As the crust thickness depends on the 
compactness of the star, we explore the allowed parameter space for the unknown high-density EOS by selecting several 
tabulated high-density EOSs that yield neutron star models with significantly different macroscopic properties. 

In the present work, we derive the perturb ation equations for torsional oscillations of magnetized stars in general relativity 
following the approach bv lMessios et al.l pOOlr ) and neglecting spacetime perturbations (Cowling approximation). We specialize 
the magnetic field to a dipole structure, which allows us to reduce the system of equations to only one dimension. A large 
number of torsional modes (fundamental and overtones) is obtained for selected magnetar models differing in the high-density 
EOS, the crust EoS and mass. Our study is valid for moderate magnetic field strengths, up to roughly 10^'' G, assuming 
that the magnetic field is confined to the crust and neglecting the possible presence of a thin fluid ocean. The reason for 
this is twofold: on one hand, for significantly stronger magnetic fields the distortion of the equilibrium shape of the star as 
well the magnetic coupling to higher harmonics should also be taken into account; on the other hand, the Alfven velocity 
ua = B /{Aivpy^'^ « 3 X fO^ cm/s _Bi5 Pi^^^ becomes comparable to the speed of shear waves in the crust when the magnetic 
field reaches values of 10^^-lO^^G, which implies that global magnetosonic waves will be strongly coupled to shear waves and 
the torsional oscillations may no longer be confined to the crust. The effec t of including global magnetosonic waves ev en at 



smaller magnetic field strengths should also be studied in more detail, see (jGlampedakis et al.ll2006l : ISotani et al 



2006l l. 



The article is structured as follows: in the next two sections we describe the general-relativistic ideal MHD equations, the 
background stellar configuration and the way that we have derived the dipole geometry of the magnetic field. In the fourth 
section we derive the perturbation equations for torsional oscillations of magnetized stars in the Cowling approximation, 
while in the fifth section we present the obtained frequencies for selected magnetar models along with fitting formulae for the 
infiuence of the magnetic field in the oscillation spectrum. The article closes with a summary and discussion. 

Unless otherwise noted, we adopt units of c = G = 1, where c and G denote the speed of light and the gravitational 
constant, respectively, while the metric signature is ( — , -|-, +). 



2 GENERAL-RELATIVISTIC IDEAL MHD 

The energy momentum tensor T'^" for a magnetized relativistic star in equilibrium, is the sum of the stress-energy tensors of 
a perfect fiuid T""'"-''' and the magnetic field T""^'^'' 

jiM" _ rpt^i^ipf) _|_ J,^^u{M) 

where 

yM^P/) ^ _^pgA^-^ (9) 

rp^-^{M) ^ H^^u" + ^H''g'"' ~~ H^H^ (10) 

Above, e, p, and H'^ are the energy density, the pressure, and a normalized magnetic field, which we define by absorbing 
a factor of \/47r, i.e., H'^ = B'^/\/A^, and = H^H^. Here and throughout the paper we will assume the ideal MHD 
approximation. In addition, we assume the shear modulus to be isotropic, so that there is no contribution of the shear stress 
in equilibrium. 

Using the stress-energy tensor (|8j , one obtains the equations of motion of the fluid by projecting the conservation of the 
energy-momentum tensor on to the hypersurface normal to it'', i.e., 

h^^T% = Q, (11) 
where = g'^^ + u^u^ is the projection tensor. The equations of motion are 

(6 + P + i/2)KV = -/i'"^ (p+i//^) +h''^{H-^H-).^^. (12) 

Due to the ideal MHD approximation the electric field 4-vector vanishes, iJ^ = F^^^u" = so that the electric field is zero 

for a comoving observer. In this approximation, the Maxwell's equations Fy^^.^^i^ — 0, where Fi_ii. is the Faraday tensor, can be 
written in the simple form: 

{u^H" -u''H'').^ = 0. (13) 
In deriving the above equation the following definitions are used 

Ffi„ = u^,E^ — u^Efi ~ efivai3U°'B'^ , (14) 

B^ = ^ef^^apu" F"'^ , (15) 
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From Equation (|13|l one derives the magnetic induction equation 

H^.^u = -u'.^H'' + u%H'' + H"uc;f,u^u'', (16) 

= [a^ + uj^ - |ne) + H^u^.pu^u'', (17) 
where cr,i^ is the rate of shear tensor, tj^^ is the twist tensor and B is the expansion, defined as 

(18) 

+ U^iah^ij^) - ^©^Mf. (19) 

— Uv:ah°'f^^ . (20) 

In a following section, the above equations will be linearized, in order to obtain the perturbed equations of motion describing 
small-amplitude oscillations. But first, we describe the construction of the background equilibrium models. 



e 































3 EQUILIBRIUM CONFIGURATION 



A strongly magnetized relativistic star has a non-spherical shape, due to the tension of the magnetic field, which is expressed 
via the off-diagonal components in the stress-energy tensor T'"'. However, the deformations from spherical symmetry induced 
by the magnetic field are small for neutron stars with magnetic field usually assumed for magnetars, because the energy of 
the magnetic field is considerably smaller then the gravitational energy, i.e., 

2 



magnetic energy 



gravitational energy GNP /R 



10" 



B 



1016 [G] 



(21) 



For this reason, we neglect the deformations due to magnetic fields in the construction of equilibrium models. Since magnetars 
are also extremely slowly rotating (the intense magnetic field spins down the star on a short timescale) we can also neglect any 
rotational deformations, so that the equilibrium models can be considered as spherically symmetric solutions of the well-known 
TOV equations described by a metric of the form 



I 2 2* 1,2 , 2Aj 2 I 2/ ifl2 , . 2 

as = — e at + e ar + r [do + sm 



(22) 



where <E> and A are the function of the Schwarzschild radial coordinate r. The four-velocity of the equilibrium model is thus 

= (e~*, 0,0,0). (23) 

We supplement the equilibr ium model by a dip ole magnetic field (but see Braithwaite fc Spruit ( 20061 ) for other assump- 
tions), following the approach of lKonno et al.l 1 19991 ) i.e. we consider an axisymmetric, poloidal magnetic field, which is created 
by a 4-current — (0, 0, 0, J^). In ideal MHD the electromagnetic 4-potential yl^ is very simple and has only one component 
in spherical polar coordinates 



= {0,0,0, A4,), 
with 



(24) 



(25) 



{ Bocauet et ahll 19951 : iKonno et all 19991 ) . I'^or the metric chosen above. Maxwell's equations -F''^^ = 47rJ'' lead to the following 
elliptic equation for 



Q^2 



+ ■ 



1 d'A. 



dr 



1 cos 6 dAtf, 
sin 9 89 



= -47rj0 . 



(26) 



The 4-vectors and can be expanded in vector spherical harmonics, so that 

A^{r,9) = ae^,{r)sm9doPi,Acos9), (27) 
Mr,0) = ji^,{r)sm9doPi,Acos9), (28) 

where (r) and ^ (r) are functions of the radial coordinate only. This leads to a one- dimensional form of Equation (|26|l 



_2A d^ae^ 



+ ($'-A') 



-2A rfa^M ^MjlM + 1) , 

dr 



(29) 



For a dipole magnetic field, i.e., Im = 1, an analytic exterior solution for a\ exists, setting jj' 
1 - 2M/r (jWasserman fc Shapiroll 19831 ) 
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Figure 1. Magnetic field lines of a typical equilibrium model. At the surface, an analytic solution for an exterior dipole magnetic field 
is smoothly joined to a numerical interior solution. The grey area shows the crust region, to which the magnetic field is assumed to be 
effectively confined when computing the torsional oscillations. 



3^6 2 



8M3 



In 1 - 



2M^ 



2M 2M^ 
+ — + — 



(30) 



where fib is the magnetic dipole moment for an observer at infinity. In the interior the potential is found by solving 
numerically Equation (|29p . The current distributio n .7'i'"\ which is the s ource term of Equation (|29p must satisfy the inte- 
grability condition that the solution is stationary (iBonazzola et al. 1993h . Here we adopt the same, simplified form for the 
current distribution as the one used in lKonno et alJ 1 1999f ). i.e., 



jf)=/or^(6 + p), 

where /o is an arbitrary constant. Moreover, regularity at the center of the star implies that a^^'"' has the form 

where Oc is a constant, which is determined by the junction condition with a^^^^ at stellar surface. Finally, the components 
of the vector J/^ describing the normalized magnetic field are given by the following set of equations 

„A „A „„„ a 



(31) 



(32) 



Hr 



Anr^ sin 9 



-deA 



-ai 



-drAs = -- 



-ai. 



(33) 



(34) 



/47r sin 6 V 47r 

Figure [T] shows the magnetic field lines of the complete solution, for a typical equilibrimum configuration. We emphasize that, 
even though the magnetic field configuration has been computed throughout the star, the boundary conditions we will assume 
for the torsional oscillations will neglect the presence of the magnetic field outside the crust region, so that, effectively, that 
magnetic field will be assumed to be confined to the crust only. 



4 PERTURBATION EQUATIONS 
4.1 Covariant form 

The perturbation equations are derived by linearizing the equations of motion (|12|l and the magnetic induction equation 
(|16[) . Since we only consider torsional oscillations, which are of axial type and do not induce density variations in spherical 
stars, there is no significant variation in the radiative part of the metric describing the gravitational field. In addition, 
since torsional modes are essentially material velocity oscillations, the imaginary part of their frequency due to the emission 
of current-multipole gravitational radiation is much smaller than the real part. For these reasons the frequency of torsional 
oscillations is determined with satisfactory accuracy even w hen neglectin g entirely the metric perturbations by setting 5g^v = 
(i.e. adopting the relativistic Cowling approximation), see LeinsI (1994) and we adopt this approximation here. 



The linearized form of the equation of motion (|12p is 
(e + p + iy^)(5<,u" = -{5e + 5p+2H^5H'^)u%u" -{e + p + H'^)u''.Ju'' 
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+ h''^ [H°'SH'' + 5H°'H'' - (7°" {5p + HpSH")] ^ - h^^ST"";,^'^ . (35) 

where ST'^"^"^ is the hnearized shear stress tensor. The latter is assumed to be related to the linearized shear tensor SS^,j, 
through 

= -2^l6S^^, (36) 
see ISchumaker fc Thornel (|l983l ). In addition, linearization of the magnetic induction equation p6|) yields 



where we used u°'Sua-f3 = —Su^Ua-p. Notice that the above relation differs in some terms from Equation (12) in lMessios et al 



(2001), due to a sign error in the latter referencqj. In Ap pendix [S] we a l so giv e the corrected equations and numerical results 
for the special magnetic field configuration considered in lMessios et al.l 1 200 ll ). 



4.2 Axial Perturbations 



The usual method for studying matter perturbations of spherically symmetric backgrounds is to decompose the perturbed 
quantities into scalar and vector spherical harmonics of definite indices I and m. In this way, the different oscillation modes 
characterized by I and m decouple and one can study them independently. Here we restrict attention to axial-type perturba- 
tions, for which the perturbed matter quantities can be written as 



5e = 5p= 0, 
5u — Su~ = 5u^ — 0, 



Su"" = 
where 



m 



sin 



'dty{t,r)b{e), 



dePi{cose), 



(38) 
(39) 
(40) 

(41) 



where de denotes the partial derivative with respect to 6. Above, y{t,r) describes the radial dependence of the angular 
displacement of the stellar material, Pi{cos9) is the Legendre polynomial of order £ and we have set m — 0, due to the 
degeneracy in m for spherically symmetric backgrounds. The above assumptions lead to the following form of the (^-component 
of the linearized equation of motion psp 



(e + p + iy^)e"*<5i 



5ff 



+ 



5Hr.4, 



2 cot OH" +[^^r + -]H' 



+ 5H'' 



H'^„ + 3cot ejy* 



8,4, 



ST 



r(p(s) 



ST' 



e</>(s) 



+ $ ^ + A r Ur 



nr0(s) 



3 cot eST 



(42) 



where we used the relation 



i^Ha, which is derived from the Maxwell's equations. 



For the derivation of the lin earized shear tensor SS^v, we adopt the relationship (5(t^^ = CuSS^v IjCarter &: Quintana 
I972I : ISchumaker fc Thornel Hiiil ) , which becomes 



Sa^^ = e SSy.v,t 

In the Cowling approximation the only nonvanishing components of (Ja^j^ are written as 

Scjr^, ~ e~'^y,trhsai? 6 and Sae^ 



1 2 1 • 2fl 

-^r e y^tb.e sm 0, 



It follows that the only non- vanishing components of the linearized shear stress tensor are 



ST 



r4, = ^A"' ^,r&sin f 
^T^B6 = ~fJ,r'yb,g sin" ( 



(43) 

(44) 

(45) 
(46) 



^ A sign error adopted initially from another publication propagated throughout the magnetic terms in the equations and final numerical 
results presented in that paper. We have confirmed that all equations and results in the non-magnetized limit were correct, while also 
the qualitative behaviour of the fundamental torsional mode as a function of the magnetic field strength was correct (see Appendix (All . 
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The linearized induction equation (|37|) yields that following relations for the first time-derivative of the components of SH'^ 

5H\t = e-''H'^5u^,t, (47) 

5H\t = -e"^ H\^5u'*' , (48) 

5H\t = -e^H^Su-*", (49) 

SH'^^t = e* [{^.rH'' - H'f'J) 5u^ + H^'Su^^^ + H''Su^,e\ ■ (50) 

Substituting the previous relations into the linearized equations of motion (|42|l and assuming that all perturbed variables 
have a harmonic time dependence, such that y(t, r) = e'"'3^(r), one obtains 



cj e y = e fiy + 



e 3^ 



r 



H';^H'^,r + H^^H"^ + A' + ^ ) + + 3 cot eH'^^H" 



2 cot + ($' + ^) + i/^ {H\,^ + + i/^^ J + H'-JJ^ 



2 cot eH^'H" +U' + i\ {H^ + H'' {H\,. - ff^^) + H'^H^g y' 



+ [H'-fy" + 



$' + - ) h^'h" + if ® , + 2 cot e(iy'')' + H'^ - 



+ 2H^H0y'Y-f + {H'Yy^-f, 



(51) 



where a prime (') denotes the derivative with respect to r. The above eigenvalue equation for the mode-frequency uj is still 
written for a general equilibrium magnetic field H^(r,6,(j)). 



4.3 Dipole magnetic field 

Restricting attention to an axisymmetric poloidal magnetic field, the eigenvalue equation (|51[) becomes 



[, + p + H'-Hr + H'Hsye-^^y = e'^y + [(^ + - A') ^ + 



+ 



2 cot en^'H" + ( $' + ± ) [H^'f + H''H\, + H^H\g y' + {H^'fy 



4 + (Hy 



y 



+ I + ^) H'-h" + H^H",,. - cot e{Hy + H'h' 
where we used the following equation satisfied by the angular function b{6) 
b,ee + 3 cot eb,e + (^ + 2)(^ - l)fe = . 



y + 2H''H"y' 



(52) 



(53) 



The above eigenvalue equation can be simplified further by assuming a dipole magnetic field of the form (|33|l . (|34p . We achieve 
this by first writing Equation (I52|l in the form 

[Ae{r) + Beir)]deYi,o + [C^r) - Be{r)] sin' OdgYeo + Vt{r) sing cos 0y,o = 0, (54) 

where Yeo is the spherical harmonic function with m = 0, which is related to P^(cos^) via the normalization Yio — 
y/{2l + l)/{ATv)Pi. The four radial functions Ae(r), Be{r), Ct{r), and Vi(r) are 

'4 



At{r) 



My + 



$' - A' U + 



y 



. ^ ^ 2 -2* (^ + 2)(£~1) 



2A-,, 



^^(0 = ^ + [(*' - «i + 2»i'] y + [a;'e-'*+'^ai + ($' - A') a/ + a/'] 3^} , 

C,(r) = ^{2aiy+ [c.W-(^ + 2)(£-l)]ai'3^}, 
0.(r) = ^^-tila,{2a/y+[($'-A')ai' + ai"]3^}, 

in which the magnetic field now appears through the radial function ai of Equations (|33|) and 



(55) 
(56) 
(57) 
(58) 
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Multiplying Equation (|54|) with ObY/q and integrating over the sphere, the eigenvalue equations takes the 1-dimensional 
form 

£{£ + 1) (A + Be) + Cf^ {Ce - Be) + C^^Vi = . (59) 
In deriving this equation we used the following properties of spherical harmonic functions: 

n^) {deYe'^) dfl, (60) 



£{l+l)Ae ~ ^Ae> I {deY^ 



Ct^Ai ~ 



-£{£ + 3)Qe+2Qe+iAe+2 + [i^Ql+i + {£ + l)'Qf] Ae - {£ - 2){e + l)QeQe-iAt-2, (61) 
^Ai' (deYlo) sine cos 9Ye,od.n 

£Qi+iQe+2Ai+2 + kQl+i -{1+ 1)Q?] A - (^ + l)QeQe~iAe-2, (62) 



\l (2^-l)(2£ + l)- ^^^^ 



From Equations (|59p . (|6ip and (|62ll it is evident that the presence of the dipole magnetic field causes the eigenfunction of 
a specific torsional {£, 0) mode to be no longer described by an angular part that corresponds to a single value of £ (as in 
the nonmagnetized case). Instead, the eigenfunction of a specific mode also acquires pieces with £ + 2, £ — 2 and higher order 
dependence (through successive couplings). In order to solve the full problem to a desired accuracy, one has to solve a linear 
system of equations of the type (|59|l for required values of £, up to some finite ^max. Since the coupling is purely a magnetic 
field effect, in the absence of a magnetic field the radial functions Be, Ce, and De vanish and the different £ dependencies 
decouple. 



4.4 Neglecting £±2 couplings 

As mentioned in the introduction, we limit our current study to moderate magnetic field strengths of a few times 10^^ G, since 
for larger values one would also have to take into account global magnetosonic oscillations. We further simplify the numerical 
problem, by neglecting the couplings to ^ ± 2 terms in Equation (I59p . This is motivated by the expectation that higher-order 
contributions will be less significant than the lowest-order piece of an eigenfunction, since the higher-order pieces are sampling 
the background magnetic field at smaller scales. This expectation may be valid for the study of low-order modes. However, it 
is difficult to predict the effect of ^ ± 2 couplings on high-order modes without solving the full two-dimensional problem. Our 
numerical results will thus be valid only to the extend that magnetic-field-induced couplings to £±2 terms can be neglected. 
Neglecting I ±2 couplings leads to the eigenvalue equation 



ai 



M + (l + 2Ai)^ 



y" + {(l + $'-A')M + M' + (l + 2A0^[($'-A')a. + 2a.']}y 



2 \ \ 
, /I , r,\ \ "1 1 2A MO-l 

e + p+ l + 2Ai — e - ^r-^ 



where 

A = l{l+l). 



Ai 



tQe+i - (^ + 1)Q? 



-(A -2) 



1(1+1) 



lie 



Aiai 
27rr'' 



■)+(2 + 5Ar)^{(<I.'-A')ai' + a/'}|^ = 0, 



{2£- l){2l + 3)' 

In order to solve a system of first-order ODEs, we define new variables yi and 3^2, through 

yi = yr'-'. 



y2 



/i + (l + 2Ai) — 



e y r 



(64) 

(65) 
(66) 

(67) 
(68) 



In defining the new variables, we took into account the form of the eigenfunction y near the center, where y 
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Torsional Oscillations of Relativistic Stars 9 




Figure 2. Mass-radius relationship of nonrotating equilibrium models constructed with various high-density EoS, in combination with 
the two choices for the crust EoS. Individual models for which we compute torsional modes are shown with symbols. 



final first-order system of equations to be solved numerically for the real eigenvalues lj is then 



3^2' 



£- 1 



yi + 



Tvr-^fi+ (1 + 2Ai)ai' 



-#-l-A-,, 

rc y2, 



(69) 



2 /2 
, /, , ON \"1 \ -2A oil \ 2 2{A-*) 

e-|-p+ (1 + 2Ai) — 7 - Aie - — - ] u re ' ' 



-(A -2) 



5Ar)^(^^-2.,.) 



e 3^1 



e + 2 



y2, 



(70) 



where Equation (|29|l was used in order to eliminate the term of ai" . 



4.5 Boundary conditions 

If one would solve the above system of equations from the center of the star to the surface, then regularity at the center 
implies the boundary condition 



y2 = ie- 1) 



^t + (l+2Ai)- 



e yi . 



(71) 



If the equations (|69p - (|70p are used for the study of global oscillations, i.e. magnetosonic waves in the core and the crust and 
shear waves in the crust, an appropriate boundary condition at the base of the fluid/crust interface will be ^'2 ' ~ 3'2^'' which 
leads to the condition 



y 



'(-) 



1 + 



y 



'(+) 



(72) 



This is analogous to the one used by Glampedakis et al. 1 20061 ) in their study of a corresponding toy problem. 

However, here, we only consider torsional modes confined to the crust and impose a zero traction condition at the base 
of the crust, which implies that 3^2=0 there. At the stellar surface, the "zero-torque-at-surface" condition is imposed, i.e. 
grj.{s)r ^ Q jSchumaker fc Thornel [liii ) . which again implies 3^2 = with our choice of variables. We neglect the possible 
presence of a thin fluid ocean. 
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Table 1. Mass (M), Radius (ij) and relative crust thickness (Ar/R) for specific equilibrium models 
constructed with various high-density EoS (A, WFF3, APR and L) in combination with two crust EoS 
(DH and NV). The subscript in the model names is indicative of the mass of each model. 



Model 


M/Mq 


i?(km) 


Ar/R(%) 


— 


Model 


M/Mq 


K(km) 


Ar/Ri%) 


A+DH14 


1.4 


9.49 


5.15 




A+NV14 


1.4 


9.48 


6.79 


A+DH16 


1.6 


8.95 


3.54 




A+NV16 


1.6 


8.94 


4.69 


WFF3+DH14 


1.4 


10.82 


6.45 




WFF3+NV14 


1.4 


10.82 


8.53 


WFF3+DH16 


1.6 


10.61 


4.98 




WFF3+NV16 


1.6 


10.61 


6.54 


WFF3+DH18 


1.8 


10.03 


3.54 




WFF3+NV18 


1.8 


10.03 


4.64 


APR+DH14 


1.4 


12.10 


7.72 




APR+NV14 


1.4 


11.93 


9.98 


APR+DH16 


1.6 


12.09 


6.23 




APR+NV16 


1.6 


11.95 


8.08 


APR+DH18 


1.8 


12.03 


5.05 




APR+NV18 


1.8 


11.92 


6.54 


APR+DH20 


2.0 


11.91 


4.05 




APR+NV20 


2.0 


11.82 


5.26 


APR+DH22 


2.2 


11.65 


3.15 




APR+NV22 


2.2 


11.59 


4.11 


L-I-DH14 


1.4 


14.66 


10.16 




L+NV14 


1.4 


13.58 


12.18 


L+DH16 


1.6 


14.78 


8.49 




L+NV16 


1.6 


13.82 


10.22 


L+DH18 


1.8 


14.83 


7.11 




L+NV18 


1.8 


14.00 


8.65 


L+DH20 


2.0 


13.82 


5.99 




L+NV20 


2.0 


14.09 


7.32 


L+DH22 


2.2 


14.73 


5.02 




L+NV22 


2.2 


14.11 


6.18 


L+DH24 


2.4 


14.54 


4.16 




L+NV24 


2.4 


14.02 


5.15 


L+DH26 


2.6 


14.13 


3.32 




L+NV26 


2.6 


13.68 


4.10 



5 NUMERICAL RESULTS 

5.1 Description of the EoS and equilibrium models 



We have studied torsional modes for a variety of neutron st ar models using four dif ferent equations of state for th e core, 
ranging from a very soft EoS (Eo S A) llPandharipande 1971) to a very stiff (EoS L ) I Pandharipande fc Smith 1975 ). with 
two intermediate ones, EoS WFF3 jWiringa et al.lll988l ') and APR (jAkmal et al.lll998l '). For each EoS we constructed a number 
of models, starting from a gravitational mass of IAMq and reaching close to the maximum mass limit in increments of 0.2Mq. 
In order to separately investigate the effect of the composition of the cru st, we matched the v arious high-density EoS to two 
different propose d equations of state for the crust, one recent derived by Douchin fc Haensell 1 200lh (DH) and, for reference, 
and older EoS by Negele fc Vautherin ( 19731 ) (NV). Figure [2] displays the mass-radius relationship of nonrotating equilibrium 
models constructed with the various high-density EoS, in combination with our two choices for the crust EoS. The individual 
models for which we compute torsional modes are shown with symbols, while their detailed properties are listed in Table [T] 

From Figure [2] it is evident that for the soft equations of state (A and WFF3) the choice of the crust EoS does not 
affect significantly the bulk properties of the star, while for the stiffest EoS L models with the same mass have considerably 
different radii. The reason for this is that the models constructed with the stiff EoS have central densities that are much 
closer to the density of the base of the crust than for soft EoS. The two crust EoS differ signifi cantly both in the detailed 
composition, as well as in the dens ity at the base of the crust , which is at p ^ 2.4 X lO^^gr/cm^ for lNegele fc VautherinI (|l973l ) 
and at p « 1.28 x 10^*gr/cm'^ for Douchin fc Haensell ( 2001 ). For the stiff EoS these different properties of the crust have a 
considerable effect on the bulk properties. 

For the two crust EoS the shear m odulus p can be computed using Equation ^ . For the NV EoS a simple fitting formula 
was already derived by Duncan ( 19981 ) 



p = 1.267 X 10^° erg cm'^ptf , 



which we also adopt here. For the DH EoS, we derive the following fitting (see also ( Pirol 2005h ) 

4/3/^\V616\4/3/l-X„x*/3 
V53J V A J 



p = 1.74 X 10''" erg cm 



-a 4/3 

Pl4 



0.41 



^ r.30 

10 erg cm 



(0.02123 + 0.37631pi4 + 3.13044pt4 



4.71841p?4 + 2.46792pi4) 



(73) 

(74) 
(75) 



where X„ is the fraction of nucleons in the neutron gas outside nuclei. 



5.2 Torsional-mode frequencies in the non- magnetized limit 

Even though the excitation of these modes and the SGR activity are likely related to the strong magnetic field, the torsional 
mode frequencies are not affected by the magnetic field for B < 10^^ G (at least when the latter is confined to the crust). 
It is not clear whether magnetars have magnetic fields that only approach this value or that are even stronger. One should 
therefore first consider torsional mode frequencies in the non-magnetized limit. 
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In Table[2]we list the frequencies of the fundamental and torsional modes tto for ^ = 2 to ^ = 10, for the equilibrium models 
discussed earlier and shown in Table [1] In addition to those models, we computed a few more i.e. L+DH17, L+NV19, and 
L+NV25 in an attempt to fit the observed frequencies to given models. Among these models, the frequency of the fundamental 
£ = 2 mode varies from 17 to 29 Hz i.e. depending on the stellar parameters it can vary by up to 30-50%. Limiting attention 
to models with M = IAA'Iq only, the variation is from 21 to 29 Hz, which still significant. Therefore, strong constraints of 
the high-density EoS could be placed by observations of the fundamental £ = 2 torsional mode in magnetars, especially if 
the mass is deduced by other means. On the other hand, the choice of the crust EoS does not seem to affect significantly the 
frequency of the fundamental torsional mode. The typical variation observed for equal mass models with different crust EoS 
is on the order of 1-5%. The same comments apply to the behavior of higher £ modes, which are nearly equally spaced for 
£ ^ S, such that a scaling law can be derived for each model. 

The above picture is drastically altered when one considers the first overtone, eti, for w hich the obtained freque ncies are 
shown in Table|31 We observe that our numerical data are in agreement with the conclusion in [Hansen fc Cioffil ( 1980|) that the 
frequencies are practically independent of the harmonic index £. A second observation is that the variations in the frequencies 
due to different choices of both the high-density and crust EoS are significant. Notice that these observations apply for the 
higher overtones efn as well. The frequencies of the first overtones vary from 500-1200 Hz, while even for models with the 
same high-density EoS and same mass the frequency varies by up to 30% when changing the crust EoS. This behavior is quite 
useful given that there are s uggestions for an observatio n in the spectrum of SGR 1802-1-14 of a torsional mode with n = 1 
with a frequency of 626.5Hz (jWatts fc StrohmaveilEooel ') and at least one more at 1837 Hz. Comparing this to our results in 
Table [3] one could immediately exclude the soft and intermediate EoS A and WFF3, if the magnetic field is such that it does 
not affect the torsional mode frequencies. The IAMq model of EoS APR with an NV description of the crust agrees with 
the observation, as does the 2.OM0 model with EoS L and NV for the crust. EoS L also seems to be compatible with the 
observation when combined with the DH crust EoS for a mass somewhat larger than I.6M0. Finally, one can observe also 
that as the stellar mass increases the fundamental mode frequency decreases while the frequency of the 1st overtone increases. 
Note that this last observation also applies for the frequencies of the higher overtones as well (also see Table [4]). 

In the tables we have marked the frequencies which are close to the observed ones (within 2%) as bold for the source 
SGR 1806-20, underlined for the source SGR 1900-H4 and italic for the source SGR 0526-66. One can observe that some of 
the stiff models fit quite well to the observational data i.e. model L-fNV25 (where 17/0 =153 Hz) for SGR 1900+14. Still, 
model A-I-DH16 fits quite well with the observed frequencies of SGR 1900-1-14 and the £ = 11 fundamental mode (not shown 
in the above tables) is exactly at 155 Hz. On the other hand, for SGR 1806-20 the models L-fDHi7 (where 15/0 — 155 Hz) or 
L-I-NV20 (where 15/0 = 152 Hz) are in good agreement with the observational data. According to our estimates, the observed 
frequencies of 18, 29, 92.5, 150, 626.5, and 1837 Hz would correspond to 2/0, 3/0, 9/0, 15/0, efi, and respectively, although 
it is difficult to explain the observational data of 26 Hz because this data is very close to the other data of 29 Hz and it may 
be difficult to explain both frequencies by crustal torsional modes only. Additionally, model APR-I-NV14 also agrees well with 
the observational data for SGR 1806-20. In the this non-magnetized limit, the currently available observational data appear 
to exclude softer EoS, such as EoS A and WFF3, if crustal torsional modes are invoked in interpreting the data. 



5.3 Implications for crust thickness 

As mentioned in the Introduction, one expects a simple relation between the thickness of the crust, Ar/R and the ratio of 
the frequencies of the fundamental mode and higher harmonics. Equation ([3]). Using our numerical results, we can test the 
validity of the expected relation and attempt to infer the thickness of the crust from the observed frequencies in several SGRs. 
We construct the following empirical formula for the crust thickness 

-5- ~ ePn—f - ■ (,'D) 

where iPn are fitting coefficients. Using our numerical results presented in previous sections we find for the most interesting 
case oi £ = 2 and for n = and n — 1 the fitted values 2/3i — 2.23 ± 0.29 for all models with a crust described by the NV 
EoS and 2/^1 ~ 2.19 ± 0.28 for the corresponding models with a crust described by the DH EoS. These values are somewhat 
larger than the ones shown in the approximate relation ([Sjl , but this should be expected, since the Equation ([Sjl was based on 
Newtonian dynamics and heuristic arguments. 

Using the above empirical relation, important conclusions can be drawn for the stellar models listed in Table 1. The 
observed frequencies of SGR 1806-20 SGR, tfo ~ 18Hz and ffi — 626.5Hz suggest a crust ratio Ar/R ~ 0.06 which combined 
with the observed frequencies favors the model L-I-NV22. Although, this conclusion is based on a single observation, it is still 
very important because it favors neutron star models with a very stiff equation of state, implying considerably larger masses 
than the typical ones. One should still keep in mind that there is an uncertainty regarding the observed frequency at 18Hz. If 
the actual frequency for the efo mode is at 29Hz then the empirical formula (|76|) implies Ar/R ~ 0.10, which suggest stellar 
models with considerable smaller mass such as APR+NV14 again, or ones with even smaller mass. 
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Table 2. Frequencies (in Hz) of the fundamental torsional modes in the non-magnetized limit (i.e. n = 
and B = 0). Values in bold type, underline, and italic type are within 2% of observational data for SGR 
1806-20, SGR 1900-f 14, and SGR 0526-66, respectively. 



Model 


I = 


2 


^ = 3 


e = 


4 


I = 


5 


I = 


6 


I = 


7 


£ = 


8 


£ = 9 


£ = 10 


A-I-DH14 


28 


5 


45.1 


60 


5 


75 


4 


90 


1 


104 


7 


119 


2 


133.7 


148.1 


A-hDHi6 


27 


2 


43.0 


57 


7 


72 





86 





99 


9 


113 


8 


127.6 


141.3 


WFF3+DH14 


26 


3 


41.6 


55 


9 


69 


7 


83 


3 


96 


8 


110 


2 


123.5 


136.8 


V V J- i- ij -L-' -L J- 1 1) 


25 


2 


39.9 


53 


5 


66 


g 


79 


g 


92 


7 


105 


5 


118.4 


131.1 


WFF3+DH18 


24 


3 


38.4 


51 


5 


64 


3 


76 


8 


89 


3 


101 


6 


114.0 


126.2 


APR-I-DH14 


24 


6 


38.9 


52 


2 


65 


1 


77 


8 


90 


4 


102 


9 


115.4 


127.8 


APR-I-DH16 


23 


4 


37.0 


49 


6 


61 


9 


73 


9 


85 


9 


97 


8 


109.6 


121.5 


APR-I-DH18 


22 


3 


35.2 


47 


3 


58 


9 


70 


4 


81 


9 


93 


2 


104.5 


115.7 


APR-hDHao 


21 


3 


33.6 


45 


1 


56 


3 


67 


3 


78 


1 


89 





99.7 


110.5 


APR-hDH22 


20 


2 


31.9 


42 


8 


53 


4 


63 


8 


74 


1 


84 


i 


94.6 


104.8 


L-I-DH14 


21 


6 


34.1 


45 


7 


60 





68 


1 


79 


1 


90 


1 


101.0 


111.9 


L-HDH16 


20 


6 


32.5 


43 


7 


54 


5 


65 


1 


75 


6 


86 


1 


96.5 


106.9 


L-HDH17 


20 


1 


31.8 


42 


6 


53 


2 


63 


5 


73 


8 


84 





94.2 


104.4 


L-HDH18 


19 


7 


31.1 


41 


7 


52 





62 


2 


72 


2 


82 


3 


92.2 


102.2 


L-HDH20 


18 


9 


29.9 


40 


2 


50 


1 


59 


9 


69 


6 


79 


2 


88.8 


98.4 


L-HDH22 


18 


2 


28.8 


38 


6 


48 


1 


57 


5 


66 


8 


76 


1 


85.3 


94.5 


L-I-DH24 


17 


5 


27.7 


37 


2 


46 


4 


55 


5 


64 


5 


73 


4 


82.3 


91.2 


L-hDH26 


17 





26.8 


36 





44 


8 


53 


6 


62 


3 


70 


9 


79.5 


88.1 


A-I-NV14 


28 


7 


45.4 


60 


9 


76 





90 


8 


105 


4 


120 


2 


134.7 


149.2 


A-hNVi6 


27 


4 


43.3 


58 


1 


72 


4 


86 


6 


100 


6 


114 


5 


128.4 


142.2 


WFF3-fNVi4 


26 


7 


42.2 


56 


6 


70 


6 


84 


A_ 


98 


1 


111 


6 


125.2 


138.6 


WFF3+NV16 


25 


4 


40.2 


53 


9 


67 


2 


80 


4 


93 


4 


106 


3 


119.2 


132.0 


WFF3+NV18 


24 


4 


38.6 


51 


7 


64 


5 


77 


1 


89 


6 


102 





114.4 


126.7 


APR-I-NV14 


25 


2 


39.8 


53 


4 


66 


6 


79 


5 


92 


4 


105 


2 


118.0 


130.7 


APR-I-NV16 


23 


8 


37.6 


50 


5 


63 





75 


3 


87 


5 


99 


6 


111.7 


123.7 


APR-I-NV18 


22 


6 


35.7 


47 


9 


59 


8 


71 


4 


83 





94 


5 


106.0 


117.4 


APR-hNV2o 


21 


4 


33.9 


45 


5 


56 


7 


67 


8 


78 


8 


89 


7 


100.5 


111.4 


APR-hNV22 


20 


3 


32.1 


43 


1 


53 


8 


64 


3 


74 


7 


85 





95.3 


105.6 


L-I-NV14 


23 


2 


36.6 


49 


2 


61 


3 


73 


3 


85 


I 


96 


9 


108.6 


120.3 


L-HNV16 


21 


8 


34.5 


46 


3 


57 


7 


69 





80 


2 


91 


3 


102.4 


113.4 


L-HNV18 


20 


7 


32.7 


43 


9 


54 


7 


65 


4 


76 





86 


5 


97.0 


107.4 


L-hNVi9 


20 


2 


31.9 


42 


8 


53 


3 


63 


8 


74 


1 


84 


3 


94.6 


104.8 


L-hNV20 


19 


7 


31.1 


41 


7 


52 


1 


62 


2 


72 


3 


82 


3 


92.3 


102.2 


L-hNV22 


18 


8 


29.7 


39 


8 


49 


7 


59 


4 


69 





78 


6 


88.1 


97.6 


L-hNV24 


18 





28.4 


38 


1 


47 


5 


56 


8 


65 


9 


75 


1 


84.2 


93.3 


L-HNV25 


17 


6 


27.8 


37 


3 


46 


5 


55 


5 


64 


5 


73 


5 


82.4 


91.2 


L-HNV26 


17 


2 


27.2 


36 


5 


45 


5 


54 


± 


63 


2 


71 


9 


80.6 


89.3 



5.4 Effect of the magnetic field 



We have already discussed earlier tlie possible effect of tiie magnetic field on tiie frequencies of the various torsional modes, 
when the magnetic field is confined to the crust . An approximate formula for the effe ct of the magnetic field based on 
Newtonian physics has been suggested by Duncan ( 1998h . see also a detailed discussion by Messios et al. ( 200lh . The formula 



suggests that the corrections in the frequency coming from the magnetic field scale as i.e. 



ijn 



1/2 



(77) 



where tf'n'^ is the frequency in the absence of any magnetic field and is a typical magnetic field strength at which magnetic 
field effects on torsional modes have become important, which we take here to be B^ = 4 x 10^^ G. In t he limit oi B « Bu, , 
this result agree s with earlier studies of the effect of the magnetic field on spheroidal modes derived by Carroll et al. (198a) 
and lUnno et al.l (| 19891 ) , where the effect of the magnetic field has been shown to have the behavior 



(78) 



where i&n is a coefficient depending on the parameters of the star (M, R and EoS) and it can be derived easily if the 
eigenfunctions of a specific mode are available. When one uses general relativity to calculate the frequencies of the torsional 
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Table 3. Frequencies (in Hz) of the first overtone of torsional modes in the non-magnetized limit 
(i.e. n = 1 and B = 0). Values in bold type are within about 2% of observational data for SGR 
1806-20. 



EoS 


£ = 


2 


e = 


3 


e = 


4 


£ = 5 


e = 


6 


e = 


7 


£ = 


8 


£ = 


9 


^ = 10 


A-I-DH14 


1206 


2 


1206 


8 


1207 


6 


1208.6 


1209 


7 


1211 


1 


1212 


7 


1214 


5 


1216 


4 


A-I-DH16 


1531 


3 


1531 


7 


1532 


3 


1533.0 


1533 


9 


1534 


9 


1536 





1537 


3 


1538 


7 


WFF3+DH14 


942 


4 


943 





943 


9 


945.0 


946 


3 


947 


8 


949 


6 


951 


5 


953 


7 


WFF3+DH16 


1101 


1 


1101 


6 


1102 


3 


1103.2 


1104 


2 


1105 


4 


1106 


8 


1108 


3 


1110 





WFF3+DH18 


1367 


2 


1367 


6 


1368 


1 


1368.7 


1369 


5 


1370 


4 


1371 


4 


1372 


5 


1373 


8 


APR-I-DH14 


761 


3 


762 





762 


9 


764.1 


765 


5 


767 


1 


769 





771 


1 


773 


4 


APR-I-DH16 


860 


2 


860 


8 


861 


5 


862.5 


863 


6 


864 


9 


866 


4 


868 


1 


869 


9 


APR-I-DH18 


965 


6 


966 





966 


6 


967.4 


968 


3 


969 


4 


970 


6 


972 





973 


5 


APR-I-DH20 


1083 


3 


1083 


6 


1084 


1 


1084.8 


1085 


5 


1086 


4 


1087 


3 


1088 


4 


1089 


7 


APR-I-DH22 


1238 





1238 


3 
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modes the relativistic effects such as the relativistic form of the sound speed, the redshift corrections and the EoS affect the 
weighting factor ea„ of equation (|78[) and have to be taken properly into account. 

In our numerical calculations of the torsional modes of magnetized neutron stars for the various EoS mentioned earlier 
we derived lists of frequencies for every model, for magnetic field strengths up to i? = lO'^^ G. By numerical fitting we have 
thus found the coefficients ean of the following empirical formula 



ifn 

ejn 



1 + iCtn 



1 1/2 



(79) 



The coefficients <a„ are listed in Table \S\ and the way that the various EoS (both for the core and the crust) affect 
the torsional frequencies becomes apparent. In general, it seems that models constructed with the DH equation of state are 
affected significantly more by the magnetic field than models following the NV EoS, for all harmonics. As a result, the large 
differences in the torsional mode frequencies between models constructed with the stiffest EoS L but with different crust EoS 
are diminished by magnetic field effects around B = 4 x lO'^^G and reversed for larger values of B. 

When the magnetic field strength is equal to _B = B^ — 4 x 10^^ G the frequencies have increased by up to 35% for 
the fundamental torsional mode (independently of the value of £) and up to 100% for the first overtone. The results suggest 
that for magnetic field strengths exceeding roughly 10^^ G the shift in the frequencies is significant and should be taken into 
account in any attempt to fit the observational data with specific stellar models. 
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Table 4. Frequencies (in Hz) of the overtone of torsional modes for 1 = 2 and n = 2, 3, and 4 in the 
non-magnetized Hmit (i.e. B = 0). Values in bold type are within about 2% of observational data 
for SGR 1806-20. 



EoS 


n = 2 


71 = 3 


71 = 4 


A-I-DH14 


2011.4 


2774.8 


3476.4 


A-I-DH16 


2552.6 


3521.8 


4412.2 


WFF3+DH14 


1572.0 


2164.7 


2708.7 


WFF3+DH16 


1835.8 


2530.7 


3167.6 


WFF3+DH18 


2278.2 


3142.1 


3933.3 


APR-I-DH14 


1270.0 


1747.9 


2186.8 


A 1~*T1 1 T~\TT 

APR-l-DHie 


1434.7 


1976.0 


2472.5 


APR-I-DH18 


1610.0 


2218.8 


2776.6 


APR-I-DH20 


1805.4 


2489.4 


3115.8 


APR-fDH22 


2062.9 


2844.2 


3559.7 


L-I-DH14 


885.1 


1217.6 


1525.1 


L-I-DH16 


978.6 


1347.9 


1688.6 


L-I-DH17 


1029.7 


1418.7 


1777.3 


L-l-DHis 


1081.5 


1490.4 


1867.3 


L-I-DH20 


1188.4 


1639.2 


2054.1 


L-I-DH22 


1314.4 


1813.3 


2272.4 


L-I-DH24 


1457.6 


2011.5 


2521.1 


L-I-DH26 


1658.5 


2289.4 


2869.5 


A-I-NV14 


1687.9 


2390.3 


2905.7 


A-I-NV16 


2113.5 


2997.8 


3675.8 


WFF3+NV14 


1315.1 


1860.7 


2250.9 


WFF3+NV16 


1536.7 


2176.5 


2647.9 


WFF3+NV18 


1898.3 


2691.5 


3295.1 


APR-I-NV14 


1094.5 


1547.3 


1864.0 


APR-I-NV16 
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1730.5 


2098.8 


APR-I-NV18 
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APR-I-NV20 
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2159.5 


2641.2 


APR-hNV22 
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3008.3 


L-I-NV14 
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1453.1 


L-HNV16 
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1316.9 


1586.8 


L-HNV18 
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1426.4 


1727.1 


L-I-NV19 


1047.1 


1482.6 


1799.9 


L-HNV20 


1091.6 


1546.0 


1879.5 


L-HNV22 


1184.1 


1678.3 


2048.6 


L-HNV24 


1294.1 


1835.3 


2248.6 


L-hNV25 


1366.2 


1937.6 


2375.4 


L-hNV26 


1462.3 


2074.3 


2546.7 



In Figures[3]and|4]we show some examples of the effect of the magnetic field on the torsional mode frequencies. Notice that 
for B < Bfj, the effect of the magnetic field on the frequencies follows a quadratic increase, in agreement with the approximate 
relation (|78p . On the other hand, when B > the modes change character and become dominated by the magnetic field, 
while the frequencies tend to become less sensitive to the stellar parameters. In this regime, additional effects, such as the 
coupling to global magnetosonic waves and to higher-order harmonics should be taken into account. 



6 SUMMARY AND DISCUSSION 

We have derived the formalism for computing torsional oscillations of relativistic stars endowed with a strong dipole magnetic 
field, confined to the crust. Our equations are valid in the Cowling approximation (no spacetime perturbations) and we 
have neglected the effect of the magnetic field on the equilibrium configuration and on the coupling of torsional modes to 
£ ±2 terms and to global magnetosonic modes. Under these approximations, our formalism allows us to obtain an estimate 
of the magnetic field effects on torsional modes, up to moderate values of the magnetic field strength (up to a few times 
Bf^ = 4 X 10^^ G). These moderate values of the magnetic field strength are appropriate for models of magnetars, for which 
there is strong evidence that they are at the heart of the SGR phenomenon. 

We have done a systematic search of parameter space by computing torsional mode frequencies for various values of the 
harmonic index £ and for various overtones, using an extended sample of models of compact stars. These models vary in 
mass, high-density equation of state and crust model, uniformly covering the allowed mass vs. radius parameter space. Our 
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Table 5. The values for the fitting factors ia„ of equation 
l|79|l . The fitting factors have been calculated for magnetic 
field strength up to 10^^ G. 



Model 


2^0 


2"! 




501 


8^0 


8"! 




0.38 


1.44 


0.67 


1.63 


0.70 


1.65 


A+DHie 


0.39 


1.34 


0.70 


1.52 


0.73 


1.55 


WFF3+DHi4 


0.45 


1.55 


0.70 


1.75 


0.72 


1.77 


WFF3+DHi6 


0.42 


1.43 


0.68 


1.62 


0.71 


1.64 


APR+DH14 


0.54 


1.63 


0.71 


1.84 


0.72 


1.87 


APR+DH16 


0.50 


1.53 


0.70 


1.73 


0.72 


1.75 


APR+DH18 


0.48 


1.44 


0.70 


1.63 


0.72 


1.66 


APR+DH20 


0.48 


1.37 


0.70 


1.55 


0.72 


1.57 


APR+DH22 


0.47 


1.32 


0.72 


1.49 


0.75 


1.51 


L+DH14 


0.80 


1.83 


0.77 


2.06 


0.75 


2.09 


L+DH16 


0.73 


1.68 


0.73 


1.90 


0.73 


1.93 


L+DH18 


0.67 


1.59 


0.72 


1.80 


0.72 


1.83 


L+DH20 


0.63 


1.50 


0.71 


1.70 


0.71 


1.73 


L+DH22 


0.60 


1.44 


0.71 


1.64 


0.72 


1.66 


L+DH24 


0.57 


1.38 


0.71 


1.56 


0.72 


1.59 


L+DH26 


0.55 


1.32 


0.72 


1.50 


0.74 


1.53 




0.26 


0.56 


0.32 


0.63 


0.33 


0.64 


rt.-riN V It) 


0.24 


0.50 


0.32 


0.57 


0.33 


0.58 


WFF3+NV14 


0.34 


0.60 


0.33 


0.68 


0.33 


0.69 


WFF3+NV16 


0.30 


0.55 


0.33 


0.62 


0.33 


0.63 


APR+NV14 


0.42 


0.63 


0.35 


0.72 


0.33 


0.73 


APR+NV16 


0.38 


0.58 


0.34 


0.66 


0.33 


0.67 


APR+NV18 


0.35 


0.55 


0.33 


0.62 


0.33 


0.63 


APR+NV20 
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0.52 


0.33 


0.59 


0.33 


0.60 


APR+NV22 


0.30 


0.49 


0.34 


0.56 


0.34 


0.57 


L+NV14 
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0.37 


0.78 


0.34 


0.79 


L+NV16 


0.51 
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0.36 


0.73 


0.33 


0.73 


L+NViB 


0.47 


0.60 


0.35 
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0.33 


0.69 


L+NV20 


0.44 


0.57 


0.34 


0.64 


0.33 


0.65 


L+NV22 
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0.54 


0.34 


0.61 


0.33 


0.62 


L+NV24 


0.39 


0.51 


0.34 


0.58 


0.33 


0.59 


L+NV25 


0.38 


0.50 


0.34 


0.57 


0.34 


0.58 


L+NV26 


0.36 


0.49 


0.35 


0.56 


0.34 


0.57 



numerical results have shown that torsional mode frequencies are sensitive to the crust model if the high-density equation of 
state is very stiff (such as EoS L). In addition, torsional mode frequencies are drastically affected by a dipole magnetic field, 
if the latter has a strength exceeding roughly 10^^ G. The effect of the magnetic field is surprisingly sensitive to the adopted 
crust model. Using our extended numerical results we have derived empirical relations for the effect of the magnetic field on 
torsional modes as well as for the crust thickness. We compare our numerical results to observed frequencies in SGRs and 
find that certain high-density EoS and mass values are favored over others in the non-magnetized limit. On the other hand, if 
the magnetic field is strong, then its effect has to be taken into account in attempts to formulate a theory of asteroseismology 
for magnet ars. Thi s topic, as well as the inclusion of global magnetosonic modes will be discussed in a separate publication 
( Sotani et ahlbood ). 
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Figure 3. The frequencies of the fundamental n = and the first overtone for f = 2, 3 and 4 torsional modes as functions of the 
normalized magnetic field (B/Bfj,). The neutron star mass is IAMq and we show results for EoS A+DH. 
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APPENDIX A: MONOPOLE MAGNETIC FIELD 

If the magnetic field is assumed to be axisymmetric. Maxwell's equations (|13|l become 
H\, + H\e + (^A' + ^) + cot oh" = 0. 

As a toy model, we consider a monopolar magnetic field in the crust 
= 0. 



For this case, the solution of Maxwell's equation (|A1|I we is (|Messios et al.ll200H ) 



where C is some constant. For this case the perturbation equation (I52p is rewritten as 

(^ + 2)(£-l) 2A 



t I I ur rr \ 2 — 2* + 2A-,, 



+ 



{^ + WHr)y" 

4 



- + - A' U + m' + - MWHr 
r ' 



y 



Ate 3^ 



If we introduce the new variables, such as 

= y, 

y2 = {fi + H^Hr)e^-''yi', 



then equations HA5|) reduce 

y.' 



fl + H^Hr 



e-^+^y^, 



y2' 



(^ + 2)(^- 1) 2* , ^ ^„r„. 2 



e yi — y2 

r 



(Al) 

(A2) 
(A3) 

(A4) 

(A5) 

(A6) 
(A7) 

(A8) 
(A9) 



In Figure \KT\ we show the fundamental and first overtone frequencies for neutron star models with mass 1.4M0. It is 
apparent that in the case of a monopole magnetic field the frequencies (especially of the fundamental mode) will be marginally 
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Figure Al. The frequencies of the 2^0 and 2^1 torsional modes as functions of the normalized magnetic field (B/Bf^ ) for EoS A+DH 
and A+NV are shown. The magnetic field has monopole geometry while the neutron stars have masses IAMq. 



affected by the magnetic field and only for stengths on th e order of 10 G its influence becomes si gnificant. This is in q ualitative 
agreement with the conclusions in lMessios et al.l ( 200 Ih . In addition, we confirm the finding in lMessios et al.l 1 200lh that the 
frequency of the fundamental mode decreases for a monopole magnetic field, as is apparent in Figure lAll 



